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Abstract 

Background: CYP19A1 encodes the aromatase which catalyzes the final reaction of estrogen biosynthesis. The 
bovine genome also contains a non-coding copy of CYP19A1, the transcribed pseudogene CYP19P1. Whereas 
CYP19A1 is transcribed in all estrogen-producing tissues, mainly in the placenta and gonads, the CYP19P1 transcript 
so far was detected in the placenta. Strikingly, one sequence segment of both transcripts exhibits an exceptional 
high identity of 98%, which implies selective pressure and suggests some kind of function. Only recently, indeed, 
coding-independent functions of several transcribed pseudogenes were reported. Therefore, we analyzed CYP19P1 
and CYP19A1 transcripts with the aim to detect clues for gene-pseudogene interference. 

Findings: The CYP19P1 transcript was first examined in silico for the presence of microRNA coding sequences 
and microRNA targets. Further, to identify tissues where CYP19P1 and CYP19A1 transcripts are co-expressed, as a 
pre-requisite for transcript interference, expression profiling was performed in a variety of bovine tissues. Our in 
silico analyses did neither reveal potential microRNA coding sequences, nor microRNA targets. Co-expression of 
the CYP19 loci was demonstrated in placental cotyledons and granulosa cells of dominant follicles. However, in granulosa 
cells of dominant follicles the concentration of CYP19P1 mRNA was very low compared to CYP19A1 mRNA. 

Conclusions: CYP19P1 and CYP19A1 transcripts might interfere in placental cotyledons. However, in granulosa cells of 
dominant follicles relevant interference between gene and pseudogene transcripts is unlikely to occur because of the 
very low CYP19P1/CYP19A1 transcript ratio. 
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Findings 

Background 

CYP19A1 encodes the aromatase which catalyzes the 
conversion of androgens to estrogens. During an earlier 
screening of a bovine placental cDNA library for CYP19A1 
clones, we also isolated cDNA clones of a homologous 
pseudogene, CYP19P1 [1,2]. Gene and pseudogene tran- 
scripts exhibit major differences due to the loss of several 
exons in the CYP19P1 transcript. Furthermore, numerous 
mutations in the pseudogene sequence produced multiple 
translational stop codons in all reading frames and thereby 
abrogated its protein-coding function. Strikingly, however, 
mutations are unevenly distributed. A sequence segment 
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of 177 bp corresponding to exon 5 of CYP19A1 is highly 
conserved, showing 98% identity, compared to 89% se- 
quence identity in general [2]. Both CYP19 loci are located 
on the same strand of chromosome 10, being separated by 
20 kb of genomic DNA [3,4]. 

Pseudogenes for long were considered as defunct copies 
of functional genes. However, recent evidence suggests 
that some pseudogenes might exert coding-independent 
functions [5,6]. Interestingly, gene-pseudogene interfer- 
ence was demonstrated by selective knock-down of the 
ABCC6P1 pseudogene, which led to a decreased transcrip- 
tion of its ABCC6 parent gene [7]. This recent evi- 
dence prompted us to take up again our analysis of the 
CYP19P1 pseudogene with the aim to detect clues of 
gene-pseudogene interference. To this end, transcripts 
were searched in silico for microRNA- coding sequences 
and microRNA targets. Further, expression profiles of 
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CYP19A1 and CYP19P1 were analyzed in a variety of 
bovine tissues. 

Materials and methods 

Samples from placentas, ovarian granulosa cells, fetal 
ovaries, endometria, adrenal glands and livers were 
collected from slaughtered cows in a local abattoir. 
Tissue samples were stored in RNAlater (Qiagen, Hilden, 
Germany) at -20°C. Granulosa cells were frozen in liquid 
nitrogen and stored at -80°C. Placental cotyledons and 
caruncles were separated manually. Despite careful separ- 
ation, caruncle samples might contain traces of coty- 
ledonary cells. Ovarian dominant and pre-ovulatory 
follicles collected before and after the LH-surge, re- 
spectively, were identified as described in [8]. Follicles 
were punctured with 18G needles and granulosa cells were 
aspirated. Total RNA was prepared using the NucleoSpin 
RNA II Kit (Macherey-Nagel, Duren, Germany), according 



to the supplier's protocol. This procedure included on- 
column DNasel digestion to remove traces of DNA. RNA 
was quantified in a NanoDrop 1000 spectrophotometer 
(PeQLab, Erlangen, Germany). RNA integrity was con- 
firmed by denaturing agarose gel electrophoresis. 

The abundance of CYP19P1 and CYP19A1 transcripts 
was measured by quantitative reverse transcription PCR 
(qPCR). For normalization purposes, the RPLP0 tran- 
script encoding a ribosomal large subunit protein was 
also measured. RNA (100 ng) was reverse transcribed 
in a 25 ul reaction volume using a mixture of random 
hexameric and oligo dT primers (Roche, Mannheim, 
Germany) and M-MLV reverse transcriptase (Promega, 
Mannheim, Germany). cDNA was purified with the High 
Pure PCR Product Purification Kit (Roche). For subsequent 
real time PCR, cDNA was amplified in a 12 ul reaction 
volume with the SensiFast SYBR No-ROX Kit (Bioline, 
Luckenwalde, Germany) using following primer pairs: 
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Figure 1 Expression of CYP19P1 and CYP19A1 in bovine tissues. A) Quantitative reverse transcription PCR analysis of CYP19P1 and CYP19A1 
transcripts. Concentration values were normalized using RPLP0 as an internal control. The abundance of the RPLP0 transcript in the bovine tissues 
is indicated in the inserted diagram. The results are shown as means +/- SEM of n = 3 independent experiments. Analyzed tissues were placental 
cotyledons (Cot) and caruncles (Car), fetal ovaries (fOv) and granulosa cells from dominant and pre-ovulatory follicles (dGC and pGC, respectively). 
CYP19A1 and CYP19P1 transcripts were not detected in pGC, adrenal glands, endometria and livers. B) Agarose gels confirming the generation of 
only the expected products during the qPCR analysis shown in the panel A. The products of qPCR reactions including reverse transcription are 
indicated below the gels by white boxes labeled + RT. Control PCR reactions without prior reverse transcription (indicated by black boxes labeled -RT) 
did not yield products. Lanes labeled with +, - and M contain positive PCR controls (PCR products from cloned cDNAs), negative PCR controls 
(PCR reactions without templates) and molecular weight markers, respectively. 
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CYP19Pl_iov, 5 '-TCATTACAACGCATCCCCAGGTT 
GA-37CrP19Pl_rev, 5 '-CTAGGTCCATGACGGGCT 
GGTATCA-3', CYP19Al_ior, 5 '-GGATCGGCAGTG 
CCTGCAATTACTA-37CyP19i41_rev, 5'-ATGCCGA 
TGAACTGCAACCCAAGTT-3 ' and RPLP0_for, 5'- 
TGGTTACCCAACCGTCGCATCTGTA-37RPLP0_rev, 
5 '-CACAAAGGCAGATGGATCAGCCAAG-3 ' (Sigma- 
Aldrich, Taufldrchen, Germany). The primer pairs were 
designed to flank intronic sequences of genomic DNA. 
The expected products from cDNA were 156, 174 and 
140 bp in length, respectively. The amplification and 
quantification of resulting PCR products were performed 
in a Light- Cycler 480 instrument (Roche) under following 
cycling conditions: Pre-incubation at 95°C for 5 min, 
followed by 40 cycles of denaturation at 95°C for 20 s, 
annealing at 60°C for 15 s, and extension at 72°C for 



15 s, and single point fluorescence acquisition at 75°C 
for 10 s to avoid quantifying primer artifacts. The gen- 
eration of only the expected products was confirmed by 
melting curve analysis and agarose gel electrophoresis, 
which did not reveal PCR products from genomic DNA 
(Figure IB). External standard curves were generated by 
co-amplification of various dilutions of cloned PCR prod- 
ucts (5 x 10~ 12 to 5 x 1(T 16 g DNA/reaction) with the cor- 
responding primer pairs. 

The statistical analyses were performed with the Sigma 
Plot 12.0 Analysis System (Jandel Scientific, San Raffael, 
CA, USA). 

Results and discussion 

The remarkably high conservation of CYP19P1 in the exon 
5-homologous sequence segment (henceforth referred to 
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Figure 2 Analysis of the putative CYP19P1 promoter. The upper panel depicts a part of the bovine chromosome 10 (included in the GenBank 
sequence NW_003 104282) with CYP19A1 and CYP19P1 depicted as arrow-headed boxes. The stippled line signifies 20 kb of intergenic DNA. 
Promoters (P2, P), exons (E2, E) and introns (In) are highlighted by hatching, black and white staining, respectively. Transcription start sites are 
represented by black arrows. The CYP19A1 GenBank sequence Z69242 (black horizontal lines below the chromosome 10 schematic) covers the 
proximal promoter P2, exon 2 with the translation start codon, and part of intron 2. By BLAST analysis, two homologous regions were found in the 
NW_0031 04282 sequence (shaded areas). The lower panel shows aligned sequences of homologous CYP19A1 and CYP19P1 regions. Asterisks indicate 
identical bases, dashes are arbitrarily inserted for optimal alignment. Transcribed sequences are printed in bold and transcription start sites are marked 
by arrows. The translation start codon is underlined and printed in italics. Above the CYP19A1 promoter sequence binding motifs of transcription factors 
GATA, SF1 and a TATA box are shown. 
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as Pl-exon 5) implies selective pressure. Hence, it might 
well be that the pseudogene exerts a coding-independent 
biological activity via the Pl-exon 5. Because pseudogene- 
derived small interfering RNAs were shown to regulate 
gene expression in mouse oocytes [9], we examined if also 
the CYP19P1 transcript encodes microRNAs. We per- 
formed in silico analyses of the Pl-exon 5 using the free 
software RNAfold [10] to detect putative stable pre- 
microRNA hairpins. However, no such structures were 
predicted by the program. The cellular abundance of 
the tumor suppressor gene PTEN transcript was found 
to be regulated by the transcript of the highly homolo- 
gous PTENP1 pseudogene via competition for micro- 
RNA binding [11]. To evaluate if the CYP19 gene pair 
could also interfere this way, we searched the Pl-exon 
5 for microRNA target sites using the MIRANDA software 
[12]. However, probable targets of known microRNAs were 
not found. Other pseudogenes are transcribed in an anti- 
sense orientation and lead to silencing of their parent genes 
by translational interference [5]. However, this mode of ac- 
tion can not apply to CYP19P1 and CYP19A1 which are 
both encoded by the same strand of chromosome 10 [4] 
and hence are transcribed in the same orientation. 

For possible interference, CYP19P1 and CYP19A1 tran- 
scripts need to be co-expressed. So far, expression of 
CYP19P1 has not been studied in bovine tissues other than 
placenta. To evaluate co-expression as a prerequisite for 
transcript interference, we performed qPCR to measure 
the abundance of CYP19P1 and CYP19A1 transcripts in a 
variety of bovine tissues. These included granulosa cells 
from dominant follicles and placental cotyledons, which 
were known to express CYP19A1 [13]. Further, granulosa 
cells from pre-ovulatory follicles, fetal ovaries, placental 
caruncles, endometrium, adrenal gland and liver were also 
analyzed, to find out whether the tissue-specificity of 
CYP19A1 expression is retained in CYP19P1. The results 
of the qPCR experiments are presented in Figure 1A. As 
expected, cotyledons and granulosa cells from dominant 
follicles expressed considerable amounts of the CYP19A1 
transcript. In contrast, the CYP19A1 transcript concen- 
tration was low or undetectable in the remaining tissue 
samples. Likewise, the CYP19P1 transcript was found in 
cotyledons and granulosa cells from dominant follicles, 
however, at a very low concentration. The ratio CYP19P1I 
CYP19A1 of mean transcript concentrations was 0.16 and 
0.016, respectively. Considering these results CYP19P1 
and CYP19A1 transcripts might interfere in placental cot- 
yledons, although the underlying mechanism remains un- 
clear. On the other hand it is unlikely that the CYP19P1 
transcript plays a major role in the regulation of CYP19A1 
expression in granulosa cells of dominant follicles. Inter- 
estingly, however, in granulosa cells from pre-ovulatory 
follicles CYP19P1 expression is shut down, just as known 
from CYP19A1, as a consequence of the pre-ovulatory 



LH-surge [14,8]. The apparent congruent tissue-specificity 
of CYP19 gene and pseudogene could be explained by 
promoter conservation. A genomic DNA sequence cover- 
ing the entire CYP19A1 and CYP19P1 loci is present in 
the database [GenBank:NW_003 104282]. We defined the 
genomic sequence immediately upstream of the tran- 
scribed CYP19P1 sequence as potential pseudogene pro- 
moter and compared it with promoters of CYP19A1 by 
BLAST analysis. Thereby we found, that 483 bp of the 
proximal CYP19A1 promoter 2 ([15]; [GenBank:Z69242]) 
were duplicated during pseudogene genesis (Figure 2). We 
have verified the sequence by sequencing cloned PCR 
products (data not shown). Gene and pseudogene pro- 
moters exhibit 83% sequence identity. Although binding 
motifs of GATA and SF1 transcription factors are retained, 
the putative pseudogene promoter has lost the functional 
TATA box, which might explain the low promoter activity 
and the use of an alternative transcription start site. 

Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

MC implementation, analysis of data and preparation of the manuscript. RF 
initiation of the study, analysis of data, preparation of the manuscript. Both 
authors read and approved the final manuscript. 

Acknowledgements 

We thank Dr. R. M. Brunner and F. Hadlich for providing assistance with in 
silico analyses, Maren Anders and Veronica Schreiter for excellent technical 
assistance. This work was supported by a grant from the Deutsche 
Forschungsgemeinschaft (DFG). 

Received: 30 September 2013 Accepted: 13 June 2014 
Published: 20 June 2014 

References 

1. Vanselow J, Furbass R: Novel aromatase transcripts from bovine placenta 
contain repeated sequence motifs. Gene 1995, 154:281-286. 

2. Furbass R, Vanselow J: An aromatase pseudogene is transcribed in the 
bovine placenta. Gene 1995, 154:287-292. 

3. Goldammer T, Guerin G, Brunner RM, Vanselow J, Furbass R, Schwerin M: 
Chromosomal mapping of the bovine aromatase gene (CYP19) and an 
aromatase pseudogene to chromosome 10 and syntenic group U5. 
Mom Genome 1 994, 5:822-823. 

4. Brunner RM, Goldammer T, Furbass R, Vanselow J, Schwerin M: Genomic 
organization of the bovine aromatase encoding gene and a 
homologous pseudogene as revealed by DNA fiber FISH. 
Cytogenet Cell Genet 1998, 82:37-40. 

5. Pink R, Wicks K, Caley D, Punch E, Jacobs L, Carter D: Pseudogenes: 
pseudo-functional or key regulators in health and disease? RNA 201 1, 
17:792-798. 

6. Pei B, Sisu C, Frankish A, Howald C, Habegger L, Mu X, Harte R, 
Balasubramanian S, Tanzer A, Diekhans M, Reymond A, Hubbard T, 
Harrow J, Gerstein M: The GENCODE pseudogene resource. 
Genome Biol 2012, 13:R51. 

7. Piehler AP, Helium M, Wenzel J, Kaminski E, Haug K, Kierulf P, Kaminski W: 
The human ABC transporter pseudogene family: evidence for transcription 
and gene-pseudogene interference. BMC Genomics 2008, 9:165. 

8. Nimz M, Spitschak M, Schneider F, Furbass R, Vanselow J: Down-regulation 
of genes encoding steroidogenic enzymes and hormone receptors in 
late preovulatory follicles of the cow coincides with an accumulation of 
intrafollicular steroids. Domest Anim Endocrinol 2009, 37:45-54. 

9. Tarn 0, Aravin A, Stein P, Girard A, Murchison E, Cheloufi S, Hodges E, Anger M, 
Sachidanandam R, Schultz R, Hannon G: Pseudogene-derived small interfering 
RNAs regulate gene expression in mouse oocytes. Nature 2008, 453:534-539. 



Chwalisz and Furbass BMC Research Notes 2014, 7:378 
http://www.biomedcentral.eom/1756-0500/7/378 



Page 5 of 5 



10. Gruber A, Lorenz R, Bernhard S, Neubock R, Hofacker I: The Vienna RNA 
websuite. Nucleic Acids Res 2008, 36:W70-W74. 

11. Poliseno L, Salmena L, Zhang J, Carver B, Haveman W, Pandolfi P: 

A coding-independent function of gene and pseudogene mRNAs 
regulates tumor biology. Nature 2010, 465:1033-1040. 

12. Enright A, John B, Gaul U, Tuschl T, Sander C, Marks D: MicroRNA targets in 
Drosophila. Genome Bioi 2003, 5:R1. 

13. Vanselow J, Furbass R, Rehbock F, Klautschek G, Schwerin M: Cattle and 
sheep use different promoters to direct the expression of aromatase 
cytochrome P450 encoding gene, Cyp 19, during pregnancy. 
DomestAnim Endocrinoi 2004, 27:99-1 14. 

14. Voss AK, Fortune JE: Levels of messenger ribonucleic acid for cytochrome 
P450 17alpha-hydroxylase and P450 aromatase in preovulatory bovine 
follicles decrease after the luteinizing hormone surge. Endocrinology 1993, 
132:2239-2245. 

1 5. Furbass R, Kalbe C, Vanselow J: Tissue-specific expression of the bovine 
aromatase encoding gene uses multiple transcriptional start sites and 
alternative first exons. Endocrinology 1997, 138:2813-2819. 



doi:1 0.1 1 86/1 756-0500-7-378 

Cite this article as: Chwalisz and Furbass: Evaluation of 
coding-independent functions of the transcribed bovine aromatase 
pseudogene CYP19P1. BMC Research Notes 2014 7:378. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at f~\ RiftMM i rpntral 

www.biomedcentral.com/submit momea central 



